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Abstract 



To estimate the diffusion constant D of particles in a plasma, we develop a 
method that is based on the mean free path A for scatterings with momentum 
transfer q> T. Using this method, we estimate A and D for squarks and 
quarks during the electroweak phase transition. Assuming that Debye and 
magnetic screening lengths provide suitable infrared cutoffs, our calculations 
yield A < 18/T and D < 5/T for both squarks and quarks. Our estimate of 
A suggests that suppressions of charge transport due to decoherence of these 
strongly interacting particles during the electroweak phase transition are not 
severe and that these particles may contribute significantly to electroweak 
baryogenesis. 
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I. INTRODUCTION 



In recent years, many scenarios incorporating extensions of the Minimal Standard Model 
(MSM) have been proposed for the generation of the Baryon Asymmetry of the Universe dur- 
ing the electroweak phase transition. A much studied mechanism uses the transport of some 
quantum number, via the CP violating interactions of certain particles with the expanding 
broken phase Higgs wall, into the unbroken phase during a first order electroweak phase 
transition. This transported quantum number, such as the axial top quark number, then 
biases the equilibrium in the direction of baryon number violation, ultimately generating a 
baryon asymmetry through weak sphaleron processes |],^|]. Of the extensions of the MSM 
used in these scenarios, those based on supersymmetric models such as the Minimal Super- 
symmetric Standard Model (MSSM) are the best motivated. By including supersymmetric 
particles, extra CP violation can be naturally introduced into a mechanism for baryogenesis. 
It is argued that a strongly first order phase transition can be achieved if the broken phase 
right-handed stop mass is less than, or of the order of, the top quark mass In this 

case, squarks with soft supersymmetry breaking masses m s ~ T, where T ~ 100 GeV is the 
temperature of the plasma, could have an important role in the charge transport mechanism 
0. However, to study charge transport via stops, one needs to have an estimate of their 
diffusion properties in a plasma at the electroweak scale. 

In Refs. a set of approximations in conjunction with the Boltzmann equation for 

quarks in the plasma of MSM particles were used to estimate the quark diffusion constant D q . 
As strong interactions dominate the diffusion process, and since stops are strongly interacting 
particles, it has been assumed that the estimate D q ~ 6/T of Ref. || is applicable to stops as 
well (the validity of this assumption is not a priori obvious because of the different statistics, 
masses, and couplings of squarks and quarks). This estimate is derived using an approximate 
gluon propagator with the thermal mass of the longitudinal gluons m g as an infrared cutoff 
and ignoring the different thermal properties of the transverse and longitudinal gluons. The 
use of the approximate gluon propagator can only yield the leading-logarithmic behavior 
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and does not result in the correct leading a 2 s non-logarithmic contribution ||[7||- However, 
the diffusive process is expected to be dominated by the t-channel gluon exchange diagrams, 
and for these diagrams, the leading logarithm contribution is expected to be dominant ||. 
A more comprehensive treatment in Ref. J7J gives D q ~ 3/T. Again, this estimate is at the 
level of the leading logarithm. 

To study the diffusion of particles in a plasma, one must consider scattering processes 
in which the momentum transfer q is not small. For the t-channel processes we consider in 
this paper, we approximate this effect by an infrared regularization of q such that q ^ T, 
since the typical momenta of the scattering particles are of order T. That is, we assume 
that such a transfer of momentum in the scattering process randomizes the momentum 
of the diffusing particle, as an approximation of the physics invovled. This momentum 
randomization approximation is implemented naturally by using as cutoffs the longitudinal 
and transverse thermal masses of the exchanged gluons which are comparable in magnitude 
to the temperature T of the plasma. We note here that the thermal masses of gluons depend 
parametrically on the strong coupling constant g s , and that these masses are comparable to 
the temperature only for realistic values of g s ~ 1. If we take the limit in which j s < 1, 
these masses will be small compared to the temeprature T, and cannot be used as cutoffs 
in our approach. In this paper, we use the above momentum randomization approximation 
to calculate the elastic^ mean free path A associated with the diffusive processes (for which 
q > T) and relate it to the diffusion constant D. Henceforth, the words "mean free path" 
refer to this diffusive mean free path. 

In electro weak baryogenesis scenarios that use charge transport, the CP violating interac- 
tions of the charge carriers with an expanding Higgs wall eventually result in the generation 
of baryon number. However, within the Higgs wall, multiple scatterings in which the final 



* Non-elastic processes in which the species of the particle changes do contribute to the mean free 
path; however, these are not considered to be diffusive processes here. 
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momentum of the charge carrier differs significantly from its initial momentum wash out 
the asymmetry caused by CP violation and suppress baryogenesis. This effect is known as 
decoherence. These same processes also contribute to the diffusion of the particles within the 
plasma. Therefore, the diffusive mean free path A which we calculate is a relevant parameter 
for estimating the supression due to decoherence. 

In this paper, we use the method described above to estimate the mean free path A s and 
the diffusion constant D s of stops that have a soft supersymmetry breaking mass m s ~ T 
in the unbroken phase of the electroweak plasma. We use the same method to estimate the 
mean free path \ q and the diffusion constant D q for quarks and compare our values with 
those of Ref. ||. We find that our method reproduces the results of Ref. for the set 
of parameters used therein. In general, our results suggest that the values of A and D of 
squarks are close to those of quarks. 

In calculating A and D, we consider only strong interactions, for they dominate the 
diffusion of squarks and quarks in the plasma. In the case of squarks, we further assume 
that scatterings from and via the heavy gluinos and squarks do not contribute significantly, 
leaving only quarks and gluons as the dominant scatterers and mediators. As consideration 
of more scatterers can only decrease the calculated values of A and therefore D, the inclusion 
of only quark scatterers will yield an upper bound (and even a reasonable order of magnitude 
estimate) for the size of the squark diffusion constant D s , in light of the results of Ref. J?|. 
In computing D q , to facilitate comparison, we follow Ref. || and only consider t-channel 
quark-quark scattering. 

As explained above, we implement our momentum randomization approximation using 
the longitudinal and transverse gluons thermal masses m g and mt in the plasma, referred to 
as Debye and magnetic masses, respectively, as physical infrared cutoffs for the exchanged 
gluon momentum. Since these masses depend on g s , the momentum transfer q in our calcu- 
lations is parametrically of order g s T, and the use of Hard Thermal Loop (HTL) propagator 
for the gluons is valid HQ. In the limit g s 1, this corresponds to small momentum trans- 
fer. However, we note that in our approximation, we demand q ^ T. Thus, in the g s ^ 1 
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limit, in our approach, we have to abandon the thermal masses as cutoffs, and choose a 
cutoff q cut parametrically independent of g s , and such that q cut > T. Then, in this technical 
sense, in our approach, the use of the HTL gluon propagator would not be justified, since 
the momentum carried by the exchanged gluon would not be parametrically soft, that is 
0(g s T), as required by the HTL formalism Our method does not yield the leading 

logarithm behavior, as g s — > 0, obtained in Refs. ||[7| for the t-channel processes considered 
therein. Nonetheless, for the physical values of g s at the electroweak phase transition, we 
obtain similar results, suggesting that our scheme reasonably approximates the physics. For 
realistic values of g s , we use the gluon thermal masses as cutoffs in our momentum random- 
ization approximation, since they naturally arise in the plasma. Then, parametrically, we 
are allowed to use the HTL propagator for the gluon. In this work, to incorporate thermal 
effects qualitatively, we approximate the effect of the HTL propagator by separating the 
gluon propagator into transverse and longitudinal parts that in general have different ther- 
mal masses. Whereas m g is calculable at one loop, m t is not calculable perturbatively and 
is unknown. Therefore, we will present our results for two representative values of m t . 

In the next section, we describe our method for calculating the diffusion constant of 
particles in the plasma. In Section III, we present our estimates for A s , D s , X q , and D q , 
followed by a discussion of our results. The appendix contains some information on the 
approximate thermal gluon propagator we use in our calculations. 

II. CALCULATION OF THE MEAN FREE PATH AND THE DIFFUSION 

CONSTANT 

Let us consider a two body scattering process where the initial and final particles have 
4- momenta (p, k) and (p', k'), respectively. We refer to each particle by its 4-momentum for 
the rest of this section. The p-particle, whose diffusion constant we calculate, scatters from 
the /c-particle. For processes relevant to the calculation of the diffusion constant D, the final 
state p'-particle is of the same species as the initial p-particle. The p-particle, fc-particle, 
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p'-particle, and fc'-particle have thermal distributions p p , p k , p p >, and p k i, respectively. The 
density per unit volume of a particle with 4-momentum p is given by p p d 3 p/(2ir) 3 . 

The transition probability for the above process per unit volume and per unit time is 
d 3 p d 3 k d 3 p' d 3 k' 



V 



(2vr) 3 2p° (2vr) 3 2k° (2vr) 3 2p>° (2vr) 3 2k'° 



(2tt) 4 5 (4) (p+ k-p'- k')\M\' 



x p p p k (1 ±p V f) (l±p k >), 



(2.1) 



where M. is the amplitude for the scattering, and ± is for final state bosons or fermions, 
respectively. Let da be the differential cross section for this process, 



da 



(27r) 3 2p'° (27r) 3 2fc /0 
Comparing Eqs. ( |2.1| ) and Q2.2| ) yields 

d 3 p 



4 \/(P • k Y ~ m l m l 



(2.2) 



V 



(27t) 3 p° 



d 3 k 



(2ir) 3 k° 



^(p- k) 2 - m 2 p m 2 k p p p k da. 



(2-3) 



To get the rate of collision per unit time r/ 1 ) of one p-particle in the plasma, we divide 
r] by p p d 3 p/(2ii) 3 , the volume density of p-particles: 

d 3 k 



V 



(i) 



- yj \p ■ k) 2 - m 2 p m 2 k p k da. 



(2.4) 



(27r)Vfc° 

To calculate the mean free path associated with the above process, we need to find the total 
rate of collision per unit time r]^J t {p) for one particle with initial 4-momentum p into any 
final state in the allowed phase space, using the total cross section a. From Eqs. ( |2.2|) and 
(fZ4|) we get 

2 r d 3 k r d 3 p'd 3 k' 



Vtot(p) 



- S W(p + k-p'- k')\M\ 2 p k {l ± Ppl ){\ ± p y ). (2.5) 



(4tt) 5 J p°k°J p'°k'° 
Note that each scatterer included will give an additive contribution of this form to rjl^[. 

The collision time r(p), the length of time between two successive collisions for a p- 
particle, is the inverse of 77^ (p) 

1 



r{p) 



Vtotip) 
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(2.6) 



The distance l(p) such a p-particle travels between two collisions is then given by l(p) = 
(\p\/p°) t(p). We finally get the mean free path A for the p-particle by taking the thermal 
average of l(p), using the thermal distribution of the p-particles. We thus get 



A 



d 3 p 



d 3 p \p\ 



(2.7) 



where we have used Eq. (|2.6|) . Note that this mean free path vanishes if the cross section 
suffers from infrared divergences. However, for diffusive processes, these divergences are 
suitably regulated and only processes with nontrivial momentum transfer contribute. The 
resulting mean free path Q2.7|) can then be related to the diffusion constant D by the relation 



D = - Xv, 
3 ' 



(2.8) 



where v is the mean velocity of the diffusing particle. 



III. RESULTS AND DISCUSSION 

We begin this section by describing some of the thermal properties of gluons in the 
plasma and how we incorporate these properties into our calculations. Due to interactions 
with the plasma, gluons develop temperature dependent masses. The longitudinal gluons 
have a thermal Debye mass m g (T) = ^8na s T, where a s = gl/ (An), at the 1-loop level and 
the transverse gluons have a non-perturbative thermal magnetic mass m t (T) that is zero 
at the 1-loop level and is expected to be O(g^T). Thus, we may assume that the infrared 
screening of longitudinal gluons occurs at a momentum scale m g , and the similar scale for the 
transverse gluons is likely lower. At the electroweak phase transition temperature T c pa 100 
GeV, a s ~ 0.1 and m g pa 1.6 T. Since the magnetic mass is unknown, the choice of transverse 
infrared momentum cutoff is rather arbitrary. However, because g s ~ 1 at scale T c , it is 
reasonable to assume that mt is of order T. In Table [I], we take as two representative values 
m t = T and m t = m g . 

The amplitude for the t-channel squark-quark diagram of Fig. [I] is 
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FIG. 1. i-channel squark-quark Feynman diagram. 
M sq = -8tt a s T a jfD^qik^T^qik), (3.1) 

where T a is a generator in the adjoint representation of the SU(3) C color gauge group, and 
q is a quark spinor. We work in the Landau gauge where, as explained in the appendix, 
[k' — k)^D^ u = 0. Since p+p' = 2p+k — k', M. sq does not depend on j/ explicitly. We use the 
amplitude M. sq of Eq. (|3.1| ) to estimate the squark mean free path A s and diffusion constant 
D s from Eqs. ( |2.7| ) and (|2.8|). The amplitude is squared and summed over all 72 species 
of quark scatterers (3 colors, 2 spins, 6 flavors, and antiparticles). We have numerically 
computed \ s and D s for supersymmetry breaking squark masses ranging from 50 GeV to 
200 GeV and observed that their mass dependence is weak. Our results are presented in 
Table |, where the values of X s and D s have been computed for m s = 100 GeV. 

In Ref. §, D q is computed using only the t-channel quark-quark scattering amplitude 
of Fig. |2|, given by 

M qg = -A7ia s q(p')T a Yq(p)D^q(k')T a Yq(k). (3.2) 

To compare our method with that of Ref. 0, we have computed D q , again numerically, 
using only M. qq (and again summing over 72 species of scatterers). The entries in Table | 
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FIG. 2. i-channel quark-quark Feynman diagram. 

labeled " JPT" refer to the numbers we get for the set of parameters that are used in Ref . || , 
namely a s = 1/7 and m g = m t = ^8na s T = 1.9 T. We see that using the JPT parameters, 
we obtain the estimate D q ~ 6/T of Ref. ||, where only the leading logarithmic contributions 
were considered. This suggests that our momentum randomization approximation, where the 
momentum transfer q > T, yields reasonable estimates for the diffusion constant. However, 
we also note that in Ref. || the diffusion constant D ~ l/ffi'f ln^" 1 )], where the logarithmic 
dependence comes from the use of a g s dependent m g as the infrared regulator. In our 
approach, as mentioned before, for physical values of g s , we can use m g , and perhaps m t , as 
the infrared cutoffs, but in general, the cutoff we choose is not a function of g s , and does not 
vanish as g s — > 0, for it has to be chosen to satisfy q ^ T. In this way, our computations 
will yield D ~ 1/ (gf) , which represents the a s dependence of the amplitudes fl37[|) and ( |3.2| ) . 
Therefore, the behavior of our results differs from those of Ref. by 1/ ln^" 1 ). However, 
for physical values of g s , the logarithm is of order unity. Taking the result of Ref. as 
a fair estimate of D q , we believe that our results are reliable up to factors of order unity. 
Note that we do not consider all the processes that contribute at this level: scatterings from 
on-shell gluons in the plasma also provide a substantial contribution. However, the results 
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5/r 
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14/T 


4/T 


JPT 



TABLE I. Results for A and D 



of Ref. |7[ suggest that the inclusion of these diagrams will not change the results by more 
than a factor of 2. 

Our calculations suggest that \ s ~ \ q and D s ~ .D 9 , up to factors of order unity, and 
most likely to within 30%, and that the effects of different statistics, masses, and couplings 
on the values of D and A for squarks and quarks are not strong. For mt = T < m g , we 
roughly get A < 10/T and D < 3/T. On the other hand, if m t = m g , our results increase 
by about a factor of 2. 

Electroweak baryogenesis scenarios that use quarks or squarks for charge transport in a 
first order phase transition suffer from a suppression due to decoherence that is caused by 
multiple scatterings of these strongly interacting particles across the width of the expanding 
Higgs wall. A measure of the strength of decoherence is the ratio of the mean free path of 



the particles to the width w of the wall 0. In Ref. [pLQl , a 2-loop MSSM calculation of the 
Higgs wall profile g 1VGS IV ~ 25/T. Using this value of w, our results yield X/w ^ 1/2, which 
suggests that although the effects of decoherence are not negligible, they are not severe. 
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APPENDIX: THERMAL GLUON PROPAGATOR 

In this appendix, we give the expression we use for the approximate thermal gluon 
propagator in a plasma, taking the different properties of the longitudinal and transverse 
gluons into account as represented by their respective cutoffs m g and m t . Let = (1, 0, 0, 0) 
be the 4-velocity of the plasma in the plasma frame. We denote the 4-momentum q of 
the propagating gluon by q^ = (q°, q) in the plasma frame. The component of n that is 
orthogonal to q is given by h, where 

n" = n" - {n ' q ^ ■ (Al) 
q l 

We define two projection operators P^ v and Q^ v , where 

p - n n ^ nu (\?) 

r fxv — <Jfiu ^2 ~ 2 { ^ - 

and 



= (A3) 

The expression for the Landau gauge thermal gluon propagator in our approximation is 
then given by 



A (L) = 



±P +-Q 



(A4) 



where q\ = q 2 — m 2 and ql = q 2 — m 2 and from q^P^ = q^Q^u = 0, we have q^D^J = 0. 
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